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Abstract. We present results from the spectral analysis of the Seyfert 1 .9 galaxy MCG -5-23-16, based on AS CA, BeppoSAX, 
Chandra and XMM-Newton observations. The spectrum of this object shows a complex iron Ka emission line, which is 
best modeled by a superposition of a narrow and a broad (possibly relativistic) iron line, together with a Compton reflection 
component. Comparing results from all (six) available observations, we do not find any significant variation in the flux of both 
line components. The moderate flux continuum variability (about 25% difference between the brightest and faintest states), 
however, does not permit us to infer much about the location of the line-emitting material. The amount of Compton reflection is 
lower than expected from the total iron line EW, implying either an iron overabundance or that one of the two line components 
(most likely the narrow one) originates in Compton-thin matter. 
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1. Introduction 

X-ray spectra of type 1 Active Galactic Nuclei (AGNs) have 
been extensively studied by many satellites in the last 15 years. 
After Ginga it has become clear that typical X-ray broad band 
spectra of Seyfert galaxies result from the superposition of dif- 
ferent components. The most common features arising over the 
primary power law emission are two reprocessed components: 
the so-called Compton reflection component at high energies 
(above ~10 keV) and the Fe Ko- emission line at 6.4 keV from 
neutral iron. Both are due to reflection of the primary radia- 
tion by optically thick, neutral or low ionized, circumnuclear 
matter. At low energies, typically below about 1 keV, a further 
component may arise (soft excess), the origin of which is still 
rather unclear 

The study of the reprocessed components and in particular 
of the iron line profile provides important information about 
the region from which they originate. If the line is produced 
in the innermost regions of the accretion disc, its profile must 
be broad, asymmet ric and double-peak ed due to kinematic and 
relativistic effects ( Fabian et al. 2000, and references therein), 
while a narrow profile indicates an origin far from the nu- 
cleus, either in th e torus envisaged in the Unification model 
jAntonuccil 1 993l) or in the Broad Line Regions. 

While ASCA observations indicate tha t a relativis tic line is 
a common feature in Seyfert Is (Nandra et alJll997l) . XMM- 
Newton have so far found unambiguous evide nce for rel- 
ativistic lines in an handful of objects only dWilms et alJ 



fO?, 'Turner et al.'2002'), while finding routinely a narrow line 
ee ves 2002: Bianchi et al. 2003a). 

It is therefore very important for understanding the proper- 
ties of the innermost regions of the accretion disc to search for 
more relativistic lines, to assess their frequency and strength. 
The fact that the relativistic line is broad, and therefore not 
easy to be separated from the continuum, and the confusing 
presence of the narrow component, implies that, unless the line 
equivalent width (EW) is very large, it should necessarily be 
searched for in bright sources. 

MCG -5-23-16 is a nearby (z = 0.0083) X-ray bright nar- 
row emission line galaxy. Its 2-10 keV flux has varied by a 
factor of 4 wi thin about 10 y ears: from ~ 8 x 10 " erg cm"^ 
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1978 (Tennant 1983), to ~ 2 x 10"" erg cm" 
in 1989 ( Nandra & Pounds 19941) . and again to a high stat e 
of ~ 9 X 10"" erg cm'^ s^^ in 1996 dWeaver et alJll998h . 
The nucleus is obscured by neutral matter with a column den- 
sity of ~ 10^^ cm"^; no signature of a warm absorber is de- 
tected. A soft excess is present in the spectrum below about 1 
keV. RossiXTE detected a Compton reflection signature in this 
Seyfert galaxy for the first time (Weaver et al. 1998 ): ASCA 
showed a strong (EW~ 300 eV) and broad Fe Ka line with a 
possible complex profile ("Weaver etal."l99'A: BeppoSAX con- 
firmed the presence of the Compton reflection component and 
measured an equivalent width of ~ 100 eV for the ir on line 
when fitted with a single gaussian ("Perola et al."2002V More 
recently Chandra HETGS spectra, briefly described in Weaveil 
( I2OOII) . showed a narrow (FWHM <4000 km s ') hne with 
EW~ 90 eV, while XMM-A^ewfon revealed a complex Fe Ka 
line constituted by a narrow unresolved component (EW~ 40 



2 



I. Balestra, S. Bianchi, & G. Matt: The reprocessing features in tlie X-ray spectrum of the NELG MCG -5-23-16 



Table 1. The log of all analysed observations, and exposure 
times. Values in brackets show net exposure times after filtering 
for particle induced flares. 



Date 


Mission 


Instr. 




11/05/1994 


ASCA 


SISO-1 


34 


29/11/1996 


ASCA 


sisO-1 


27 






LEGS 


36 


24/04/1998 


BeppoSAX 


MEGS 


77 






PDS 


33 


14/11/2000 


Chandra 


ACIS-S HETG 


76 


13/05/2001 


XMM-Newton 


EPIG-PN 


38 (6.0) 


01/12/2001 


XMM-Newton 


EPIG-PN 


25 (19.1) 



eV) and a broad (FWHM ~ 40 OOP km s~' ) cornponen t with an 
equivalent width of ~ 120 eV dPewangan et alJ2003h . 

In this paper we present a reanalysis of the ASCA, 
BeppoSAX, Chandra and XMM-A^ewfo« observations based 
on a model which includes Compton reflection and two dis- 
tinct components, a narrow and a broad one, for the iron line. 
We will exploit the different capabilities of the various satellites 
(energy resolution for Chandra, sensitivity for XMM-Newton, 
hard X-ray coverage for BeppoSAX) to obtain informations on 
the location and physical properties of the emitting regions. 

2. Observations and data reduction 

2.1. ASCA 

The source was observed by ASCA twice, on May 1994 and 
on November 1996 (Tab le H). Data re l ative t o the first obser- 
vation were published bv lWeaver et al . Vl997^, those relative to 
the second one were briefly mentioned by Weaver et al. ( 1998). 
The spectra were extracted from the screened event files taken 
from the Goddard archive. In this paper we will only deal with 
the SISO-1 spectra (0.5 - 10 keV). Data were analysed with 
HEAsoftS.l andXsPEC 11.2.0. 

2.2. BeppoSAX 

BeppoSAX observed this source on April 1998 (Table 
Results were published by Perola et al. (2002) . Data reduction 
follow ed the standard procedure presented bv lGuainazzi et all 
(11999^. using HEAsoft 5.1. Spectra were extracted from re- 
gions of radius 8' for the LEGS and 4' for the MECS, and 
analysed with Xspec 1 1 .2.0. A normalization factor of 0.8 was 
adopted between the PDS and the MECS, appropriate for PDS 
spectra extracted with variable rise time threshold. The norma- 
lization of the LECS to the MECS was left as a free parameter 
to account for the different time coverage. 

2.3. Chandra 

The Chandra observation was performed on November 2000 
(Table with the Advanced CCD Imaging Spectrometer 
(ACIS-S) and the High-Energy Transmission Grating 
Spectrometer (HETGS) in place. The high flux of the source 



together with a default frame time of 3.2 s resulted in a 0th 
order spectrum which is strongly affected by pileup (80% 
according to WebPIMMS'). Therefore, we will only use the 
1st order co-added MEG (0.4 - 5 keV) and HEG (0.8 - 10 
keV) spectra in this paper Data were reduced with the Chandra 
Interactive Analysis of Observations software (CIAO 2.2.1), 
using the Chandra Calibration Database (CALDB 2.10). 
Grating spectra were analysed with S her pa 2.2.1. 

2.4. XMM-Newton 

XMM-Newton observed MCG -5-23-16 twice on 13 May 
2001 and 1 December 2001 for 38 ks and 25 ks respectively 
(Table^. Both observations were performed in the Full Frame 
Mode using the blocking optical Medium Filter with the EPIC 
PN and MOS simultaneously operating. We will not deal with 
the RGS spectra in this paper because they offer poor statistics. 
Data were reduced with SAS 5.3.3. Events corresponding to 
pattern were used for the PN and 0-12 for the MOS in both 
observations. The extracted background high energy (E > 10 
keV) lightcurves clearly showed particle-induced flares, which 
were filtered adopting a 3 cr threshold above the quiescent state 
background count rate. This resulted in net exposure times for 
the PN of 6.0 ks and 19.1 ks, respectively for the observations 
of May and December. In both observations the MOS count 
rates (3.5 s"' and 3.2 s"' for MOSl) are higher than the 1% 
pileup limit for this instrument (0.7 s see table 3 of the 
XMM-Newton Users' Handbook^), therefore we will not deal 
with the MOS spectra in this paper. On the other hand, the PN 
count rates (7.5 s"' and 8.0 s"') are just within the 1% pileup 
Hmit (8 s"' for this instrument). EPIC-PN spectra (0.5 - 10 
keV) and lightcurves were extracted from a radius of 40", and 
analysed with Xspec 11.2.0. 

In the following, errors are at the 90% confidence level for 
one interesting parameter (Ax^ - 2.71). 



3. Data analysis 

3. 1. Flux and spectral variability 

We analysed six X-ray observations of MCG -5-23-16, cover- 
ing a history of seven years (see Table Q). The 2-10 keV flux 
(averaged over an entire observation) of the source has varied 
by about 25% in about 1 year (see Table|3. Large short-term 
variability is also present, as best illustrated by the BeppoSAX 
observation, which is the longest available (see Fig.^. Despite 
the remarkable flux variability, however, the hardness ratio 
(4 - 10.5 keV)/(1.5 - 4 keV) is basically constant, implying 
that no spectral variability is present at these energies. Similar 
temporal behaviours are also found in the other observations. 



' http://heasarc.gsfc.nasa.gov/Tools/w3pimms.html 
" http://xmm.vilspa.esa.es/external/xmm_user_support/ 
documentation/uhb_2. 1/XMM_UHB .html 
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Fig. 1. BeppoSAX MECS hardness ratio (4- 10.5 keV)/(l .5 -4 
keV). 

Table 2. The total flux (2-10 keV) as measured adopting the 
baseline model (described further in the text). 



Observation 


Flux 2-10 keV 




in units of 10"" erg cm"- s"^ 


ASCA 1994 


9.24 ± 0.26 


ASCA 1996 


7.45 + 0.21 


BeppoSAX 1998 


9.30 ± 0.07 


Chandra 2000 


9.66 + 0.11 


XMM 05/2001 


8.08 ±0.01 


XMM 12/2001 


7.12 + 0.02 



3.2. Spectral analysis 

We will only deal with data above 2.5 keV, because we are 
mainly interested in studying the reprocessing features, i.e. the 
iron lines and the Compton reflection component. 

Our baseline model consists of an exponentially cut-off 
power law and a reflection component from a neutral Compton- 
thick slab, isotropically illuminated by the primary radiation, 
subtending the solid angle R = Q/27r and with an inclina- 
tion angle ; with respect to the line of s ight (PEXRAV model 
in XSPEC: 'Magd ziarz & Zdziarskilll995l) . We kept the abun- 
dances of the heavy elements Z fixed to the solar value. We 
also included in the model two components of the iron line 
to account for its complex profile (Figure |2j, as suggested 
both by ASCA jWeayeretal . 1997) and XMM-Newton re- 
sults jPewangan et alJ 20031) : a narrow gaussian centred at 
6.4 keV and a relativistically broadened iron line, produced 
by the accretion di sc suiTounding a Schwarzschild black hole 
jpabian et al.ll98sl) . This choice will be further justified below. 

The DISKLINE model parameterizes the radial line emis- 
sivity as a power law, i.e. oc r^i. We fixed the line rest energy 
to 6.4 keV, q - -2 and the inner radius ri„ to 6 rg (last stable 
orbit), while we kept the outer radius rout as a free parameter 
The inclination angle of the disc / has been always linked to the 
same value of the PEXRAV model. This is equivalent to assume 
that both the broad line and the reflection component are emit- 
ted in the disc. (For this reason we also tried to fit the reflec- 



Energy (keV) 

Fig. 2. EPIC-PN residuals around the iron line energy for the 
second XMM-Newton observation when fitting with a simple 
absorbed power law (upper panel), and with a power law plus a 
narrow Gaussian whose flux is fixed to 4.5 x 10"^ ph cm"^ s"', 
as explained in the text (lower panel). 



tion component with the model REFSCH instead of PEXRAV, as 
the former includes relativistic effects. Both for the BeppoSAX 
and XMM-Newton spectra, however, no significant differences 
are found in the values of the parameters and of the;i'^ between 
the two models. Therefore, and for the sake of simplicity, we 
will adopt the PEXRAV model throughout this work). 

Finally we included a uniform cold absorber with column 
density Nh to take into account local absorption, in addition to 
the Galactic one Ng = 8.12 x lO^o cm^^ (HEASARC WSnH^). 

The analysis is made a bit complicated by the fact that, 
while the XMM-A^ewfo« observation is definitely the best one 
to constrain the iron line properties, the BeppoSAX observa- 
tions is the only capable to measure accurately enough the 
reflection component. We therefore decided to adopt an iter- 
ative procedure, consisting first in estimating the inclination 
angle from the relativistic line in the XMM-Newton second, 
and longer, observation, then deriving R, with / fixed, from the 
BeppoSAX spectrum, and finally coming back to the XMM- 
Newton spectrum to check if, with the so found value of R, the 
inclination angle is still the same. Once that was verified, with 
R - 0.45 and / = 42° obtained from such procedure, we fit- 
ted all spectra with R and / fixed to these values. In performing 
the described procedure we found rout to be very loosely con- 



http://heasarc.gsfc.nasa.gov/cgi-bin/Tools/w3nh/w3nh.pl 
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Fig. 3. Chandra HEG residuals when fitting with a simple ab- 
sorbed power law around the iron line energy. 
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strained. We therefore fixed it in all fits to the value giving the 
lowest i.e. 400 rg. 

On the other hand, if ri„ is left free to vary (with rout fixed to 
400 rg), no improvement in the;^'^ is found, further justifying 
our choice of keeping rin fixed to the innermost stable orbit 
(6rg). 

The significance of the relativistic disc line is discussed in 
lOewangan et alJ ( l2003l) . Here, suffice it to say that fitting the 
spectrum with a narrow gaussian line results in x^/d-o.f. = 
243.4/158; leaving the width of the line free to vary, the fit 
improves giving ;i'^/t/.o./. = 198.4/157 (significant at more 
than 99.99% confidence level according to the F-test), with cr = 
0.28 + 0.07 keV (coiTesponding to a FWHM of 30 800 + 7 700 
km s '). Such a large width is hard to explain other than being 
due to relativistic and kinematic effects in the innermost re- 
gions of the accretion disc (see e.g. Fabian et al. 2000). Indeed, 
a better fit (x^/d.o.f. - 192.8/157) is obtained using a rela- 
tivistic line model, with all parameters fixed to the abovemen- 
tioned values but /, left free to vary as said above. 

No significant further improvement is found adding a nar- 
row line Ot'^/^^y ^ 191.5/156). However, C/iant/raHETGS, 
thanks to their excellent spectral resolution, clearly reveal the 
narrow iron line componentat En = 6.38+0.02 keV (Figure|3}, 
which is unresolved even at the gratings resolution (FWHM< 
6 500 km s"' at 99% confidence level). Therefore, given the 
Chandra detection of a narrow component, we decided to add 
it to all fits. 

To illustrate the presence of the broad line, in Figure |2]the 
EPIC-PN residuals around the iron line energy for the second 
XMM-Newton observation are shown, when fitting with a sim- 
ple absorbed power law (upper panel), and with a power law 
plus a narrow Gaussian whose flux is fixed to 4.5 x 10"^ ph 
cm~^ s"' (lower panel; see below for the choice of this value). 
The broad line appears to be slightly asymmetric, consistently 
with the significant but not dramatic improvement in the 
when fitting with a relativistic line model rather than a broad 
Gaussian. 

The detected 2-10 keV flux variations between the different 
observations of only 25% at most (Table|2} should not result in 
any measurable variability on the amount of the reflection com- 



Fig. 4. Fe Ka narrow line flux as measured adopting the base- 
line model. 
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Fig. 5. Fe Ka broad line flux as measured adopting the baseline 
model. 



ponent. This is why we assumed R to be constant and we fixed 
it to the BeppoSAX best fit value for all fits. In this respect, it 
is reassuring that the values of F so obtained are all consistent 
with one another within the error bars. The cut-off energy Ec 
has also been assumed to be constant (Ec = 110 keV). Table|3] 
shows the results obtained with the baseline model applied to 
the six observations considered. 

In Figures0]and|5lthe flux history of the narrow and broad 
line components, respectively, is shown. Both components are 
consistent with being constant within the, admittedly large, er- 
rors. In both cases this is not surprising, given the small flux 
variations in the continuum. For the narrow component, which 
is expected to originate in distant matter, this is also expected 
on theoretical grounds. To reduce the error bars on the broad 
component, we refitted all observations with the flux of the nar- 
row component fixed to its weighted mean, i.e. 4.5 x 10"^ ph 
cm"^ s ' . New results are summarized in Tableland illustrated 
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ASCA 1994 


ASCA 1996 


BeppoSAX 1998 


Chandra 2000 


XMM 05/2001 


XMM 12/2001 


Nh (10^2 cm-2) 

r 

R 


0.89+"? 

-0.35 

1 65+°'" 
0.45* 


l-21+n1i 

-0.44 

1.69 ±0.10 
0.45* 


i.57+!;» 

-0.31 

1.73 ± 0.07 
0.45+SJ?? 


1.25+"-?' 

-0.18 

1 70+°* 
0.45* 


l-94+!!i^ 

-0.40 

1.77 ±0.07 
0.45' 


1.80 ±0.23 
1.74 ± 0.04 
0.45* 


(keV) 


110* 


110* 


iiO!« 


110* 


110* 


110* 


In (lO-^phcm-2 s ') 




4 8+-'' ^ 

^•°-3.1 


7 7+43 
'•'-5.1 


8.6 ±3.5 




4 ^+0.9 


EWn (eV) 

Id (10-5ph s ') 


58 29 

12 6+'-'^ 


50i« 
9 5+" 

^■-'-5.2 


"-'-43 

< 15.4 


70 ±28 
< 17.6 




45+' 

Q Q+2.8 
^■^-2.5 


EWd (eV) 


114+87 
-69 

355.4/371 


1Q7+89 
^"'-59 

359.1/341 


< 134 
87.8/83 


< 153 
252.6/428 


114+47 

165.9/140 


199.5/157 



NOTE.-The following parameters are fixed: iron line energy E„ = 6.4 keV; disc inclination angle ; = 42° (both in PEXRAV and disc line 
model); disc line emissivity q = -2; inner and outer disc radius r^^ = 6 Tg and r„^i = 400 respectively. 
* denotes fixed parameters. 



Table 4. Fluxes and equivalent widths for the broad line component after the flux of the narrow line component has been fixed to 
the weighted mean value. Note that BeppoSAX and Chandra measured fluxes are no longer only upper limits. 





ASCA 1994 


ASCA 1996 


BeppoSAX 1998 


Chandra 2000 


XMM 05/2001 


XMM 12/2001 


Id (lO^^ph cm-2 s"') 
EWd (eV) 
;r-/dof 


13 5+''-'* 

1 97+64 
^^'-52 

356.1/372 


10 2+'-^ 

1161? 
359.1/342 


10.8!^^ 

Q7+40 

88.4/84 


10.3+^5 

911^ 
255.4/429 


9 8+^ 
101+" 
165.9/141 


9.8+^5 

114!^? 
199.6/158 



o 



A XMM 
X CHANDRA 
□ BsppoSAX 
O ASGA 



Fig. 6. Fe Ka broad line flux as measured fixing flie flux of the 
narrow line component to the weighted mean value (4.5 x 10"^ 
ph cm"^ s"'). 

in Figure |6l No significant changes are found for the other pa- 
rameters. 



4. Discussion 

We have reanalysed the XMM-Newton observations of MCG - 
5-23-16, as well as previous ASCA, BeppoSAX and Chandra 
ones, confirming t he presence of a s ignific antly broad line al- 
ready reported by 'Dewangan et al." ('2003V A fit with a rela- 
tivistic disc profile is somewhat better than that with a broad 



Gaussian. Moreover, a so broad line (cr - 0.28 keV) is hard 
to for m other than in a relativistic disc (see e.g. Fabian et al.^ 
I2OO0I) . which we therefore consider the most likely expla- 
nation. A narrow line is also present, as clearly shown by 
Chandra, as well as a Compton reflection component, meas- 
ured by BeppoSAX. Exploiting the modest long term flux vari- 
ability, we could confidently fit the XMM-A^evvfo« spectrum 
adding these two further components (with the parameter R de- 
scribing the relative amount of reflection component kept fixed 
to the value found by BeppoSAX). 

The detection of a broad, likely relativistic line in this 
source is important. Predicted on theoretical grounds, relativis- 
tic hnes w ere found to be com mon in ASCA spe ctra of AGN 
dNandra et al.. 1997). However. iLubiiiski & Zdziarski (200i 
showed that at least part of the Unes could be ascribed to a 
narrow component from distant matter, also to be expected on 
the basis of Unification Models (e.g. Ghisellini et al.. 19941) . 
Chandra and XMM-Newton have indeed confirmed the almost 
ubiquity of the narrow component, while relativistic lines have 
been unambiguotisly detected so far only in a few cases (e.g. 
IWilms et al.ll200l[ iTurner et alJl2002l) . The discovery of a rel- 
ativistic line in MCG -5-23-16 adds therefore one more case 
to a still short list. It is worth remarking that the relativistic 
lines have been observed so far in sources in which either the 
EW of the line is very large, as in the case of MCG -6-30- 
15 and NGC 3516 (.Turner et al. 2002) . or that are very bright, 
like MCG -5-23-16. The presence of the narrow component, 
in fact, could make difficult to detect the broad component in 
not very bright source s, if its EW is modest. Even if in som e 
cases, e.g. N GC 5506 dMatt et alJl200lt iBianchi et 101200331) 
NGC 7213 (iBianchi et al."2003b^. NGC 5.548 ("Pou nds et alJ 
.2003i) and NGC 4151 (.Schurch et aL.2003y) . relativistic lines 
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with significant equivalent widths are definitely absent in high 
signal-to-noise XMM-Newton spectra, the actual frequency of 
relativistic lines is still to be assessed. 

The amount of reflection component, R = 0.45, is low 
if compared with the total EW of the lines, which is about 
170 eV re.g. iMatt et alJll99lt lOeorge & F abian"l991). One 
possible explanation is iron overabundance (M att et aljll997ii) : 
alternatively, one of the two line components may originate 
in Compton-thin matter. The obvious candidate is the narrow 
component (the broad one being likely emitted by the accretion 
disc), which therefore could be associated with the BLR rather 
than the molecular torus. (It is worth noting that the column 
density of the absorber is too low to make it a good candidate 
for producing the line (e.g. Matt et al. 2003)). This hypothesis 
cannot be directly tested, because MCG -5-23-16 is a NELG 
and therefore the BLR is supposed to be obscured in the optical 
band, and therefore no BLR line width, to be compared with the 
upper limit derived from Chandra, is available. Alternatively, 
the line can be emitted in a Compton-thin torus. If indeed the 
narrow line is emitted in Compton-thin matter, no evidence for 
a Compton-thick t orus is present in this source, similar to the 
case of NGC 7213 jBianchi et aTA2003b.) . 
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